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ABSTRACT 

We describe a new and efficient technique which we call sideband or phase- 
modulation searching that allows one to detect short period binary pulsars in obser- 
vations longer than the orbital period. The orbital motion of the pulsar during long 
observations effectively modulates the phase of the pulsar signal causing sidebands to 
appear around the pulsar spin frequency and its harmonics in the Fourier transform. 
For the majority of binary radio pulsars or Low-Mass X-ray Binaries (LMXBs), large 
numbers of sidebands are present, allowing efficient searches using Fourier transforms 
of short portions of the original power spectrum. Analysis of the complex amplitudes 
and phases of the sidebands can provide enough information to solve for the Keplerian 
orbital parameters. This technique is particularly applicable to radio pulsar searches 
in globular clusters and searches for coherent X-ray pulsations from LMXBs and is 
complementary to more standard "acceleration" searches. 

Subject headings: pulsars — stars: binary — stars: neutron — methods: data analysis 

1. Introduction 

Short orbital-period binary pulsars have proven to be sensitive laboratories that allow tests 
of a variety of physical processes including theories of gravitation. Unfortunately, except for those 
systems containing the brightest pulsars, they are notoriously difficult to detect. When the ob- 
servation time T b s is greater than a small fraction of the orbital period P or b, Doppler effects due 
to the orbital motion of the pulsar cause drastic reductions in the sensitivities of pulsar searches 
(Johnston & Kulkarni 1991). 
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Numerous groups (e.g. Middleditch & Kristian 1984; Anderson et al. 1990; Wood et al. 1991) 
have developed variations of so called "acceleration" searches in order to mitigate the loss in sensi- 
tivity caused by orbital motion. These searches can almost completely recover this lost sensitivity 
when T b s < P or f,/W by taking advantage of the fact that orbital motion causes an approximately 
linear drift of the pulsar spin frequency f sp i n . Recent acceleration searches of globular clusters have 
discovered numerous binary systems, some with orbital periods as short as ~ 90min (Camilo et al. 
2000; Ransom et al. 2001; D'Amico et al. 2001). 

Acceleration searches are effective only for T Q bs ^ Porb/W, but search sensitivities improve as 
T^ 2 . In effect, acceleration searches are limited to discovering only the brightest binary pulsars in 
ultra-compact — systems with orbital periods less than a few hours — orbits. For a 2 ms pulsar in 
a 1 hour orbit with a low mass white dwarf (WD) companion, the optimal observation time Tb es t for 
an acceleration search (using eqn. 22 from Johnston & Kulkarni 1991) is only ~ 300 — 600 seconds. 
If this pulsar were located in a globular cluster, where observations of 5 — 10 hours in duration are 
common, acceleration searches would be a factor of \fTobljTbest ~ 8 times less sensitive than an 
optimal (i.e. coherent) search of a full observation. 

We have developed an efficient search technique for short orbital period binary pulsars that 
complements the use of acceleration searches. Since it requires T j, s > 1.5P or b and its sensitivity 
improves as the ratio r or b = T b s /P or b increases, it allows the search of very interesting portions of 
orbital parameter space for the first time. This technique was developed and tested over the past 
several years as part of a Ph.D. thesis (Ransom 2001). A similar but independent derivation of 
many of the properties of this technique has been recently published by Jouteux et al. (2002), based 
on an earlier mention of this work in Ransom (2000). While the basic search techniques presented 
by these papers are nearly identical, aspects concerning the determination of the orbital elements 
of a candidate pulsar discovered with this technique are presented in more detail in this work. 

2. Phase Modulation due to Orbital Motion 

During observations of a binary pulsar where T f, s is longer than P or b, the orbital motion causes 
a modulation of the phase of the pulsar signal. Differences in light travel time across the projected 
orbit advance or delay the pulse phase in a periodic fashion. 

2.1. Circular Orbits 

If the orbit is circular (where the eccentricity e and the angle of periapsis to are both zero), 
the phase delays are sinusoidal. The fundamental harmonic 4 of the resulting signal as sampled at 



4 The following analysis can be applied to any harmonic h of a pulsar's signal by substituting hr sp i n for r spin and 
the appropriate harmonic phase for <f> ap in ■ 
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a telescope can be described by 

rij = n(u) = a cos [2irr spin u + 4> spin + & orb cos(27rr orb n + 4> rb)] , (1) 

where a and <p S pin are the amplitude and phase of the pulsation, & or b and 4> or b are the amplitude 
and phase of the modulation due to the orbit, r is a Fourier frequency (i.e. r = fT f, s = T b s /P), 
and u = j/N = (tj — t )/T b s is the discrete dimensionless time during the observation of the j th 
sample of N such that < u < 1. We can write $ or b and (fi or t, in terms of the three Keplerian 
orbital parameters for circular orbits (P or b, the semi-major axis x or b = a\ sin(z)/c, and the time of 
periapsis T a ) as 

^ orb — ^"^ orb J spin — t.oto 

and 



&orb — ^ft % orb f spin — 7.375 p °/ \ (2) 

*spin (,SeCJ 



0orfc = 27T to T ° + ^, (3) 
±orb ^ 

where the units of both are radians, and the pulsar's mass function fi is defined as 

[M 2 sin(i)] 3 

;i (M 1 + M 2 ) 2 ' W 

where i is the orbital inclination and Mi and Mi are the pulsar and companion masses (in solar 
masses) respectfully. 

The Discrete Fourier Transform (DFT) of the signal in eqn. 1 when & or b = can be approxi- 
mated as (see e.g. Ransom, Eikenberry, & Middleditch 2002) 

N-l 

-2-Kiru 



iV-l i 

A(r) = V rij e~ 2nijr / N ~ N / 



n(u) e- 27nru dn (5a) 

~ Q ^ c T-4>smn, f e 2vri(r spln -r)u ^ n ^) 

2 7o 

- A „ V^^-r) ^ k(^pin-r)] 

— ^coherent c / j v u W 

ft V spin r ) 

where A co h er ent = ^ spin represents the coherent Fourier response at the spin frequency r sp i n . 
Similarly, we can approximate the DFT of our phase modulated signal from eqn. 1 as 

A(r) = f e ^r spin u+<t> spin +^ orb cos(2TTr orb u+4> orb )] e ~2mru ^ u 

2 Jo 

= A coherent I e ^or b cos(2nr orb u+<t> orb ) e 2ni(r spin -r)u du _ ^ 

Jo 

Using the Jacobi- Anger expansion (e.g. Arfken 1985) 

oo oo 

e izcos(8) = i m J m (z)e im6 = J m (z)j m ( e+ ^, (7) 
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where J m (z) is an integer order Bessel function of the first kind, we can expand the integrand in 
eqn. 6b to give 

oo „! 

A(r) = A coherent Jrn^orb) / j™^ orb u+t orb+ l) ^i(r spin -r)u dn (ga) 

oo „1 

= A coherent J m ^orb)e im ^ b+ ^ / e 2 ^ r ^ +mr °^ u du. (8b) 

m=-oo ^° 

Since the integrand in eqn. 8b is identical in form to that in eqn. 5b, the Fourier transform of a 
phase modulated signal is equivalent to the Fourier transform of a series of cosinusoids at frequen- 
cies centered on the pulsar spin frequency r sp i n , but separated from r spin by mr orb Fourier bins 
(i.e. sidebands). The response is therefore 

A(r) = A coherent V J m (d> orb ) e^^+l) j^ spin+ mr orb -r) ™ [vr(r + mr orb - r)] 
m _oo 7r(^ P m + mr orb - r) 

When r orb is an integer (i.e. the observation covers an integer number of complete orbits), the 
Fourier response shown in eqn. 9 becomes particularly simple at Fourier frequencies r sp i n + sr orb , 
where s is an integer describing the sideband in question. In this case, the summation collapses to 
a single term when m = s, giving 

A (r S pi n + Sr orb ) = ^4 C o/ierent^s( < ^orb) e ls ^ >orh+ 2 ) . (10) 

In effect, the phase modulated signal produces sidebands composed of a series of sidebands split in 
frequency from r ap ; m by sr orb Fourier bins, with amplitudes proportional to J s (Q orb ), and phases 
of s (4>orb + f ) + 4> sp in radians 5 . Figure 1 shows the Fourier response for a "typical" pulsar- WD 
system in a circular 50 min orbit during an 8 hour observation. 

While the number of sidebands implied by eqn. 9 is infinite (assuming & orb / 0), properties of 
the Bessel functions — which define the magnitude of the Fourier amplitudes of the sidebands - 
give a rather sharp cutoff at ~ 2& orb sidebands (or ~ <E> or & pairs of sidebands). The maximum value 

— 1/3 

of J s (Q orb ) also occurs near s = <& orb at a value of ~ Q.&7h<& orb , which gives the phase modulation 
response its distinctive "horned" shape when $> orb is large (see Fig. 1). The average magnitude of 

— 1/2 

the phase modulation sidebands falls off more quickly as ~ 0.611<l> or6 (see Fig. 2). 

For pulsars with narrow pulse shapes and many Fourier harmonics, phase modulation produces 
sidebands split by r orb Fourier bins centered around each spin harmonic. For spin harmonic number 
h, the phase modulation amplitude is <& rb,h = ^x orb hf sp i n . Therefore, for large & rb,h, there are 
~ h times more sidebands than around the fundamental, but the average amplitudes are smaller 
by a factor of ~ hr 1 ! 2 . 



5 Note that J- n (z) = ( — l) n J„(z), and also that J„(z) can be negative. If the sideband amplitude determined 
by the Bessel function is negative (i.e. J s (z) is negative), the measured phase of the Fourier response will differ by 
7T radians from the predicted phase since the measured amplitude is always positive. 
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2.2. Small Amplitude Limit 

When <&orb "C 1, a phase modulated signal will have only three significant peaks in its Fourier 
response. These peaks are located at r sp i n and r sp i n ± r or b and have magnitudes proportional to 
M^orb) ~ 1 - $lrb/ 4 ~ 1 and J±i(,$orb) ~ ±$orfc/2 respectively. 

Middleditch et al. (1981) performed a detailed analysis of phase modulated optical pulsations 
in the small amplitude limit from the 7.7 second X-ray pulsar 4U 1626-67 in order to solve for 
the orbital parameters of the system. The separation of the sidebands provided a measurement of 
the ~ 42 min orbital period and the magnitudes and phases of the three significant Fourier peaks 
allowed them to solve for $ or b (and therefore x or t,) and 4> or \, directly. 

Small amplitude phase modulation might also be observed from freely precessing neutron stars 
(e.g. Nelson, Finn, & Wasserman 1990; Glendenning 1995). For certain geometries, the geometric 
wobble of the pulsar beam would periodically modulate the pulse arrival times and cause the 
system to appear as a compact binary. If such rapidly precessing objects exist, they would display 
a small number (~ 1 — 3) of pairs of sidebands in a power spectrum of a long observation, since the 
maximum possible phase modulation amplitude for such a system is <3? or 5 < tt/2 radians (Nelson 
et al. 1990). 

2.3. Elliptical Orbits 

Konacki & Maciejewski (1996) derived a Fourier expansion for elliptical orbital motion with 
the useful property that the magnitudes of higher order harmonics of the orbital expansion (which 
correspond to independent modulation amplitudes <& rb,h') decrease monotonically. In the limit of 
circular orbits the expansion reduces to the single phase modulating sinusoid found in eqn. 1. 

Each harmonic of the orbital modulation generates sidebands around the spin harmonic of the 
form discussed in §2.1. The sidebands created by harmonic h' of the orbital expansion are sepa- 
rated from each other by h'r or b Fourier bins. Since <& rb,h' monotonically decreases with increasing 
orbital harmonic h', the number of sidebands generated decreases while their spacing and average 
amplitudes increase. The most important aspect of this superposition of sidebands is the fact that 
the Fourier response near a spin harmonic contains sidebands spaced by r or t, Fourier bins — from 
the fundamental of the orbital expansion — just as in the limiting case of circular orbits. 

3. Pulsars with Deterministic Periodic Phase Modulation 

Figures 3 and 4 show P sp i n , Porb, and <& rb for 54 binary pulsars from the literature, all with 
orbital periods less than 10 days. A few conclusions can be reached immediately. First, radio 
millisecond pulsars (MSPs) are the most common binary pulsar. Second, there is a relatively flat 
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distribution of P or b for radio pulsars which cuts off quite dramatically near P or b ~ 0.08 d (~ 2 h) 
- the two systems with shorter orbital periods are the recently discovered accreting X-ray MSPs 
XTE J175 1-314 (Markwardt et al. 2002) and XTE J0929-314 (Galloway et al. 2002). And third, a 
large fraction of the systems are located in globular clusters. 

The fact that large numbers of relatively compact binary pulsar systems should contain MSPs is 
not surprising in light of the standard "recycling" model for MSP creation (Verbunt 1993; Phinney 
& Kulkarni 1994). This model explains not only the preponderance of pulsars with millisecond spin 
periods, but also the fact that many MSPs have low-mass (~ 0.02 — 0.4 M Q ) companion stars in 
compact circular orbits. These recycled systems make up the vast majority of the pulsars shown 
in Figure 3 and have typical values of <E> or b of 10 2 to a few times 10 3 radians. 

The cutoff in P or b at approximately two hours is almost certainly due in some part to selection 
effects. The rapid orbital motion of these systems prevents their detection using conventional 
search techniques — including acceleration searches — for all but the brightest systems (Johnston 
&; Kulkarni 1991). In fact, the only MSPs discovered with shorter orbital periods (P or b ~ 40min 
for both XTE J1751-314 and XTE J0929-314) were discovered during transient X-ray outbursts 
when the pulsed fractions were very high (Markwardt et al. 2002; Galloway et al. 2002). Since 
selection effects keep us from detecting ultra-compact systems that are not exceptionally bright, 
it is difficult to determine how much the observed lack of such systems is due to a true lack of 
compact systems. It is certain that neutron star (NS) systems with much shorter orbital periods 
do exist, as the new accreting X-ray MSPs and other ultra-compact LMXBs such as 4U 1820—303 
(P orb ~ llmin), 4U 1850-087 (P orb ~ 21min), 4U 1627-673 (P orb ~ 42min), and 4U 1916-053 
(Port ~ 50min) (Chou et al. 2001; Liu et al. 2001) have shown. 

The large fraction of known systems inhabiting globular clusters can be explained by three 
basic facts. First, due to their known location on the sky, long targeted observations of clusters 
allow fainter systems to be detected (e.g. Anderson 1992). Second, repeated observations of certain 
clusters have allowed many very weak pulsars to be detected with the help of amplification from 
interstellar scintillation (Camilo et al. 2000; Possenti et al. 2001). Finally, large numbers of binary 
pulsars are expected to have been created in globular clusters due to dynamical interactions (e.g. 
Kulkarni Sz Anderson 1996; Davies & Hansen 1998; Rasio et al. 2000). 

Recently, Rasio et al. (2000) have conducted extensive population synthesis studies for binary 
pulsars in dense globular clusters like 47 Tucanae. Their initial results imply that exchange interac- 
tions among primordial binaries will produce large numbers of NSs in binary systems. In fact, their 
results predict that large populations of binary pulsars with low mass companions (m c ~ 0.05 M ) 
should exist in many globular clusters with orbital periods as short as ~ 15 min. 

The evidence seems to suggest that a large population of ultra-compact binary pulsars should 
exist — particularly in globular clusters — which have so far eluded detection. With long obser- 
vations, these pulsars will be detected if search algorithms can identify the distinctive sidebands 
generated by the systems. The rest of this paper discusses such an algorithm. 
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4. 



Sideband Searches 



The fact that power from an orbitally modulated sinusoid is split into multiple sidebands allows 
us to increase the signal-to-noise of the detected response by incoherently summing the sidebands. 
According to one of the addition theorems of Bessel functions 



theoretically all of the signal power can be recovered in this manner for noiseless data. Practically, 
though, complete recovery of signal power is not possible for data with a finite signal-to-noise ratio 
- especially if <& or b S> 1. Such a summation would require calculation of the ~ 2<E> or b significant 
sidebands at the precise Fourier frequencies of their peaks r sp i n ± sr or &. This computationally 
daunting task would result in very large numbers of independent search trials when $> or b, r sp i n , and 
r or fe are unknown. 

Even if we could recover all of the power using a sum of the sidebands, the significance of 
the measurement would be less than the significance of a non-modulated sinusoid with the same 
Fourier power. This loss in significance is due to the fact that noise is co-added along with signal. 
The exact loss in significance can easily be calculated since the probability for a sum of m noise 
powers to exceed some power P sum is the probability for a \ 2 distribution of 2m degrees of freedom 
to exceed 2P sum (see e.g. Groth 1975, and references therein). 

In the case where only a few sidebands are suspected to be present (i.e. <& or b < 5), a brute- 
force search using incoherent sideband summing may be computationally feasible and worthwhile. 
Unfortunately, except for binary X-ray pulsars with long spin periods and possibly certain freely 
precessing NSs (see §2.2), the number of astrophysical systems with <£ or {, < 5 is probably small 
compared to those with & or b » 1 (see §3). 



For systems with Q or b S> 1, the large number of regularly spaced sidebands allows a very 
efficient detection scheme. Since the sidebands are separated from each other by r or b Fourier bins, 
they appear as a localized periodicity in the power spectrum of the original time series which can 
be detected with a second stage of Fourier analysis. By stepping through the full length power 
spectrum and taking short Fast Fourier Transforms (FFTs) of the powers, we incoherently sum any 
sidebands present and efficiently detect their periodic nature. As an added benefit, the measured 
"frequency" of these "sideband pulsations" gives a direct and accurate measurement of the orbital 
period. 

We define the following notation to represent the two distinct stages of the Fourier analysis 
for detecting a phase modulated signal. As defined in eqn. 1, the initial time series rij contains 



oo 




(11) 
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4.1. Two-Staged Fourier Analysis 
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N points. After Fourier transforming the rij (usually using an FFT of length N), the complex 
responses at Fourier frequencies r are represented by A(r) as defined in eqn. 5a. The powers and 
phases are then simply P(r) and <f>(r) respectively. The short FFTs of the P(r) used to detect 
the sideband pulsations (i.e. the secondary power spectrum) are of length M (where M <C N and 
should closely match the extent of the sideband structure in the original power spectrum) and 
generate complex responses ^(r), powers -P2M) an d phases 4>2(r). 

In order for a pulsed signal to undergo enough modulation to produce periodic sidebands, the 
observation must contain more than one complete orbit (i.e. r or b > 1). One might initially think 
that two complete orbits would be necessary to create a Nyquist sampled series of sidebands in 
the original power spectrum. However, since the sideband pulsations are not bandwidth limited, 
periodic signals with "wavelengths" less than two Fourier bins (i.e. r or b < 2) will still appear in 
the P-2,{r). Instead of appearing at Fourier frequency r = M/r or b as for signals with r or & > 2 (see 
§4.2.2), the fundamental harmonic of the response will appear aliased around the Nyquist frequency 
( r Nyq = M/2) at r a u as = M — M/r or b. While r or b > 1 is required to begin to produce periodic 
sidebands, simulations have shown that this technique realistically requires r or & > 1.5 to create 
detectable sidebands from most binary systems during reasonable observations. 

As the number of orbits present in the data (r or &) increases, so does the spacing (in Fourier 
bins) between the sidebands. Since each sideband has a traditional sine function shape in A2(r) 
(see §2.1), the full- width at half-maximum (FWHM) of a sideband is approximately one Fourier 
bin. Therefore, an increase in spacing effectively decreases the duty-cycle of the sideband pulsations 
(i.e. ~ l/r or b). Signals with small duty-cycles are easier to detect for two reasons. First, significant 
higher harmonics begin to appear in Piir) which can be incoherently summed in order to improve 
the signal-to-noise ratio. A rule-of-thumb for the number of higher harmonics present in P2(r) is 
approximately one-half the inverse of the duty-cycle (~ r or &/2). Second, as the duty-cycle decreases, 
the Fourier amplitudes of each harmonic increase in magnitude up to a maximum of twice that of 
a sinusoid with an identical pulsed fraction (see Ransom et al. 2002, for a more detailed discussion 
of pulse duty-cycle effects). This effect makes each individual harmonic easier to detect as r or b 
increases. 



4.2. Search Considerations 

As the number of orbits present in a time series increases, the signal-to- noise ratio of a de- 
tectable pulsar spin harmonic goes approximately as ~ {r or b — 1) _/3 where 0.3 < (3 < 0.4 (see 
Figures 6 and 7). In contrast, the signal-to- noise ratio of a pulsar spin harmonic detected coher- 

— 1/2 

ently goes as T bs . This leads to a rule-of-thumb for those conducting sideband searches: "The 
longer the observation the better." When making very long observations it is important to keep the 
on-source time (i.e. the window function) as continuous and uniform as possible. Significant gaps 
in the data — especially periodic or recurring gaps of the kind common in X-ray observations - 
create side-lobes around Fourier peaks of constant-frequency pulsations that a sideband search will 
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detect. It is a good idea before a sideband search to take the power spectrum of the power spectrum 
of the known window function in order to identify spurious orbital periods that a sideband search 
might uncover. 

Since phase modulation sidebands require coherent pulsations (i.e. f sp in is constant in time) in 
order to form, it may be necessary to correct the time series (or equivalently the Fourier transform) 
such that the underlying pulsations are coherent. As an example, if the telescope's motion with 
respect to the solar system barycenter is not taken into account during long observations, inter- 
stellar pulsations — and therefore any phase modulation sidebands as well — will be "smeared" 
across numerous Fourier bins due to Doppler effects. By stretching or compressing the time series 
to account for the Earth's motion before Fourier transforming (or equivalently, by applying the 
appropriate Fourier domain matched filter after the Fourier transform, see Ransom et al. 2002) the 
smeared signal can be made coherent for the purposes of the sideband search. 

If the time-dependence of f sp i n is unknown a priori, due to accretion in an LMXB, the spin- 
down of a young pulsar, or timing (spin) noise, for example, one could in principle attempt a series 
of trial corrections and a sideband search for each trial. Such a search methodology is similar in 
computational complexity to traditional "acceleration" searches. 

While §4.1 provided basic principles for detecting compact binary pulsars using the two-staged 
Fourier analysis, in the rest of this section we provide a more detailed discussion of how a search 
might be conducted. Figure 5 shows two examples of sideband searches. The left column shows data 
from an 8 h observation of the globular cluster 47 Tucanae taken with the Parkes radio telescope at 
20 cm on 2000 November 17. The signal is the fundamental harmonic of the 2.1ms binary pulsar 
47 Tuc J {fspin = 476.05 Hz, P orb = 2.896 h, r orh = 2.762, <5> orb = 120.87 radians, and m c ~ 0.03 M ) 
with a signal-to-noise ratio of rj = a/ a n . ~ 0.0054. The right column shows data from a simulated 
8 h observation similar to the 47 Tuc observation described above. The signal is the fundamental 
harmonic of the 2.0 ms binary pulsar used in Figure 1 (f sp in = 500 Hz, P or b = 50min, r or b = 8, 
&orb = 413.44 radians, and m c = 0.2 M ) with a signal-to-noise ratio one half that of the 47 Tuc J 
observation (rj ~ 0.0027). 

Once a long observation has been obtained and prepared as discussed above, we Fourier trans- 
form the data and compute the power spectrum P(r). We then begin a sideband search which 
is composed of two distinct parts: initial detection of a binary pulsar and the determination of a 
detected pulsar's orbital elements. 

4-2.1. Initial Detection 

To maximize the signal-to- noise of a detection in P2(r), we want to match both the length and 
location of a short FFT M with the width and location of a pulsar's sideband response (~ 2r or b§ rb) 
in P(r). Because a pulsar's orbital parameters are unknown a priori, we must search a range of 
short FFT lengths M and overlap consecutive FFTs. Therefore, we must choose the number and 
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lengths of the short FFTs, the fraction of each FFT to overlap with the next, and the number of 
harmonics to sum in each FFT. These choices are influenced by the length of the observation and 
the nature of the systems that we are attempting to find. 

All currently known binary pulsars with P or b < 10 d have modulation amplitudes of & or b ^ 
6500 radians (see Figure 4), with the vast majority in the range <& or b ~ 20 — 5000 radians. This range 
also includes many currently unknown but predicted "holy grail" -type systems such as compact 
MSP-black hole systems or MSP- WD systems with orbital periods of ~ lOmin. This range of <E> or b 
values is due to the fact that the short spin periods of MSPs tend to offset the small semi-major 
axes (xorb) of their typical binary systems, while binary systems with wider orbits tend to have 
longer spin periods. Multiplying this range by twice the minimum r or b of 1.5 gives a rule-of-thumb 
range of M ~ 60 — 15000 where the high end should be increased based on expected values of r or b- 
For 8 — 12 h globular cluster searches a reasonable range of values would be 64 < M < 65536 in 
powers-of-two increments, for a total of 11 different values for M. 

The choice of the overlap fraction is more subjective. To optimize a detection in P2(?), one 
should overlap the short FFTs by at least 50% in order to avoid excessive loss of signal-to-noise 
when a series of sidebands straddles consecutive FFTs. Larger overlap percentages may result in 
higher signal-to-noise, but at the expense of increased computational costs and numbers of search 
trials (even though the trials are not completely independent when overlapping). 

In Figure 5 plots a and e, the grey shaded regions show the sections of the P(r) that were 
Fourier transformed for the sideband search. The short FFTs of length M = 1024 and M = 8192 
were centered on the known pulsar spin frequency and resulted in the P^r) displayed in plots 
b and / respectively. The short FFT lengths were chosen as the shortest powers-of-two greater 
than 2r or b<fr or b- For 47 Tuc J (plot a), the "horned" structure of the phase modulation sidebands 
described in §2.1 is easily recognizable. This shape is not apparent in plot e, and in fact, traditional 
searches would detect nothing unusual about this section of the power spectrum. 

Normalization of the P2(r) results in the usual exponential distribution with mean and standard 
deviation of one for a transform of pure noise. If one first normalizes the P(r) as shown in plots 
a and e (see Ransom et al. 2002, for several normalization techniques), the resulting P2(r) can be 
normalized simply by dividing by M. 

Plots b and / of Figure 5 show the highly significant detections of the binary pulsars in P2 (r) . 
For 47 Tuc J, only one orbital harmonic is visible in since r or b < 2, but the single-trial 

significance is approximately ~ 23 a. For the simulated 50min binary pulsar, 3 significant orbital 
harmonics are detected (predicted harmonic locations are marked with a dotted line in plot /). 
When summed, the detection has a single-trial significance of ~ 15 a. 
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4-2.2. Orbit Determination 

Once a binary pulsar has been detected, examination of the P(r) and <p(r) allow the estimation 
of the Keplerian orbital parameters and the pulsar spin period P sp i n - Brute- force matched filtering 
searches centered on the estimates are then used to refine the parameters and recover the fully 
coherent response of the pulsar. 

The simplest parameter to determine is the orbital period P or b- The spacing between sidebands 
is r or b = T b s / P or b bins, implying that the most significant peak in P 2 {r) should occur at Fourier 
frequency Mjr rb- Conversely, if we measure the Fourier frequency r meas in P 2 {r) with the most 
power, we can solve for r or b using r or b = M/r meas . By using Fourier interpolation or zero-padding 
to oversample A 2 (r) and P 2 (r), we can determine r meas to an accuracy of 

3 
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where P 2 (r meas ) is the normalized power at the peak and a is the signal "purity" — a property 
that is proportional to the root-mean-squared (RMS) dispersion of the pulsations in time about 
the centroid — which is equal to one for pulsations present throughout the data (see Middleditch 
et al. 1993; Ransom et al. 2002). Remembering that P orb = T ob s/r r b and solving gives 

meas , 3T {, S 

° rb ~ M 7raMy/6P 2 (r mea Ty 
For aliased signals, the orbital period is 

T obs (M - 

1 meas 

^orb,aliased ~ , , n/r /Fd~7 V' ^ ' 

M -KaM^QP 2 (r meas ) 

For pulsars with & rb ^ 1 observed such that r or b > 2, it is often possible to measure P or b to one 
part in 10 4 , corresponding to a ~ 1 s error for many of the pulsars shown in Figure 4. 

Once P rb is known it becomes possible to measure the projected semi-major axis x or b by 
determining the total width of the sidebands in P{r). For bright pulsars like 47 Tuc J, the number 
of Fourier bins between the "horns" (A r ~ 2<£ or f,r or {,) as shown in Figure 5a can be measured 
directly and converted into x or b using 

Xorb ^ -j 7 ^ , ( 15 ) 

J spin T orb 

where f sp i n can be estimated as the frequency midway between the phase modulation "horns" . 

For weaker pulsars, where the sideband edges are not obvious (as in Figure 5e), A r can be 
estimated using numerous short FFTs. By taking a series of short FFTs of various lengths around 
the detection region of P(r) and measuring the power, "centroid" and "purity" 6 at the frequency 



The measured values of "centroid" and "purity" are estimates of a signal's location and duration in a time series as 
determined from the derivatives of the Fourier phase and power at the peak of the signal's response (see Middleditch 
et al. 1993; Ransom et al. 2002). 
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corresponding to P or b, one can map the extent of the hidden sidebands and compute an estimate 
for x orb . 

For pulsars in circular orbits, which constitute the majority of the systems in Figures 3 and 4, 
the only remaining orbital parameter is the time of periapsis passage T a . Defining the time since 
periapsis passage as Ax- o = t D — T Q and using eqn. 3 we get 

where 4> or b can be measured from the phases of the sidebands in A(r) using the following technique. 

From eqn. 10 we know that the sidebands have phases of 4> s = s (4> or b + f ) + fispin radians. 
Unfortunately, some of the measured phases will differ by it radians from those predicted by eqn. 10 
since the corresponding sideband amplitudes (the J s (&orb)) are negative (see footnote 5). If we 
predict values for J s ($orb) using our knowledge of r or b, f sp i n , and x or b, we can "flip" (i.e. add or 
subtract ir radians to) the measured phases of the sidebands that are predicted to have negative 
amplitudes. We can then estimate 4> or b simply by subtracting the phases of neighboring sidebands 
using cj) orb = <j) s+ i - <j) s - 7r/2. 

For weak pulsars each measurement of 4> orb has a large uncertainty. Fortunately, we can make 
~ 2<i> or b measurements of 4> or b by using each pair of sidebands and then determine 4> or b and its 
uncertainty statistically. Figure 5 plots c and g show histograms of the measurement of 4> or b using 
this technique. 

It is important to mention that for weak pulsars it may be difficult to determine the location 
of the sidebands in order to measure their phases. The location of the sideband peaks can be 
calculated by measuring the phase <p m eas of the fundamental harmonic discovered in ^(r). The 
number of Fourier bins from the first bin used for the short FFT to the peak of the first sideband 
in A(r) is simply (f>meas 

r or b/2ir. Subsequent peaks are located at intervals of r or b Fourier bins. 

Once estimates have been made for the Keplerian orbital parameters, we compute a set of 
complex sideband templates (i.e. matched filters for the Fourier domain response of an orbitally 
modulated sinusoid) over the most likely ranges of the parameter values given the uncertainties 
in each. These templates are then correlated with a small region of the Fourier amplitudes A(r) 
around f sp in- When a template matches the pulsar sidebands buried in A(r), we recover essentially 
all the power in that pulsar spin harmonic (Ransom et al. 2002). Figure 5 plots d and h show the 
results of just such a matched filtering operation. The Doppler effects from the orbital motion have 
been completely removed from the data and the resulting Fourier response is that of an isolated 
pulsar. 

For binary pulsars in eccentric orbits, the techniques discussed here could in principle be 
applied to each set of sidebands from the orbital Fourier expansion (see §2.3). Such an analysis 
would be much more complicated than that described here and would require a high signal-to-noise 
detection. 
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5. Discussion 

We have described a new search technique for binary pulsars that can identify sidebands 
created by orbital modulation of a pulsar signal when T & s > P or b- Sideband or phase-modulation 
searches allow the detection of short orbital period binary pulsars that would be undetectable using 
conventional search techniques. 

Figures 6 and 7 show how the sensitivity of sideband searches compares to that of acceleration 
searches for a wide range of observation times T Q b s . Acceleration search sensitivities are near- 
optimal when T b s < 0.1P or b (Johnston & Kulkarni 1991), but degrade rapidly at longer integration 
times when the constant frequency derivative approximation breaks down. Sideband searches, on 
the other hand, require T Q b s > l.hP or b in order to approach the sensitivity of optimal duration 
acceleration searches, but for longer T a b s sensitivity improves as ~ (r or b — l)^ 13 where 0.3 < (3 < 0.4. 
For targeted searches of duration ^8—12h (e.g. typical globular cluster observations), sideband 
searches are ^2—10 times more sensitive to compact binary pulsars (P or b i$ 4/ij than optimal 
duration acceleration searches. In general, one should think of acceleration and sideband searches 
as complementary to each other — sideband searches allow the detection of ultra-compact binary 
pulsars, while acceleration searches maximize the detectability of isolated and longer-period binary 
pulsars. 

The fact that sideband searches target a different portion of orbital parameter space than 
acceleration searches and yet require significantly less computation time, provides good reason to 
include them in future pulsar searches where T j, s > 30 min. In fact, the Parkes Multibeam Pulsar 
Survey (T b s ~ 35 min, Lyne et al. 2000; Manchester et al. 2001) has included a sideband search 
in an on-going re-analysis of their survey data at the cost of only a marginal increase in computer 
time. While the probability of detecting a binary pulsar with P or b < 23 min is almost certainly 
quite low, the discovery of a single such system would provide a wealth of scientific opportunity. 

Acknowledgments We would like to thank G. Fazio for supporting our research and encouraging 
us in this work. Many of the computations for this paper were performed on equipment purchased 
with NSF grant PHY 9507695. SMR is a Tomlinson Fellow, JMC is supported by NSF Grant 
9819931 and by the National Astronomy and Ionosphere Center, and SSE is supported by a NSF 
CAREER Grant. 

REFERENCES 

Anderson, S. B. 1992, PhD thesis, California Institute of Technology 

Anderson, S. B., Gorham, P. W., Kulkarni, S. R., Prince, T. A., & Wolszczan, A. 1990, Nature, 
346, 42 

Arfken, G. 1985, Mathematical Methods for Physicists, 3 rd edition (San Diego: Academic Press) 



-14- 



Camilo, F., Lorimer, D. R., Freire, P., Lyne, A. G., k Manchester, R. N. 2000, ApJ, 535, 975 
Chou, Y., Grindlay, J. E., k Bloser, P. F. 2001, ApJ, 549, 1135 

D'Amico, N., Lyne, A. G., Manchester, R. N., Possenti, A., k Camilo, F. 2001, ApJ, 548, L171 
Davies, M. B. k Hansen, B. M. S. 1998, MNRAS, 301, 15 

Galloway, D. K., Chakrabarty, D., Morgan, E. H., k Remillard, R. A. 2002, ApJ, 576, L137 
Glendenning, N. K. 1995, ApJ, 440, 881 
Groth, E. J. 1975, ApJS, 29, 285 

Johnston, H. M. k Kulkarni, S. R. 1991, ApJ, 368, 504 

Jouteux, S., Ramachandran, R., Stappers, B. W., Jonker, P. G., k van der Klis, M. 2002, A&A, 
384, 532 

Konacki, M. k Maciejewski, A. J. 1996, Icarus, 122, 347 

Kulkarni, S. R. k Anderson, S. B. 1996, in Dynamical Evolution of Star Clusters - Confrontation 
of Theory and Observations: IAU Symposium 174, Vol. 174, 181 

Liu, Q. Z., van Paradijs, J., k van den Heuvel, E. P. J. 2001, A&A, 368, 1021 

Lyne, A. G., Camilo, F., Manchester, R. N., Bell, J. F., Kaspi, V. M., D'Amico, N., McKay, N. 
P. F., Crawford, F., Morris, D. J., Sheppard, D. C., k Stairs, I. H. 2000, MNRAS, 312, 698 

Manchester, R. N., Lyne, A. G., Camilo, F., Bell, J. F., Kaspi, V. M., D'Amico, N., McKay, 
N. P. F., Crawford, F., Stairs, I. H., Possenti, A., Kramer, M., k Sheppard, D. C. 2001, 
MNRAS, 328, 17 

Markwardt, C. B., Swank, J. H., Strohmayer, T. E., Zand, J. J. M. i., k Marshall, F. E. 2002, ApJ, 
575, L21 

Middleditch, J., Deich, W., k Kulkarni, S. 1993, in Isolated Pulsars, ed. K. A. Van Riper, R. Ep- 
stein, k C. Ho (Cambridge University Press), 372-379 

Middleditch, J. k Kristian, J. 1984, ApJ, 279, 157 

Middleditch, J., Mason, K. O., Nelson, J. E., k White, N. E. 1981, ApJ, 244, 1001 
Nelson, R. W., Finn, L. S., k Wasserman, I. 1990, ApJ, 348, 226 
Phinney, E. S. k Kulkarni, S. R. 1994, ARA&A, 32, 591 

Possenti, A., D'Amico, N., Manchester, R. N., Sarkissian, J., Lyne, A. G., k Camilo, F. 2001, 
astro-ph/0108343 



-15- 



Ransom, S. M. 2000, in ASP Conf. Ser. 202: Pulsar Astronomy - 2000 and Beyond, 43 

Ransom, S. M. 2001, PhD thesis, Harvard University 

Ransom, S. M., Eikenberry, S. S., & Middleditch, J. 2002, AJ, 124, 1788 

Ransom, S. M., Greenhill, L. J., Herrnstein, J. R., Manchester, R. N., Camilo, F., Eikenberry, 
S. S., & Lyne, A. G. 2001, ApJ, 546, L25 

Rasio, F. A., Pfahl, E. D., & Rappaport, S. 2000, ApJ, 532, L47 

Verbunt, F. 1993, ARA&A, 31, 93 

Wood, K. S., Hertz, P., Norris, J. P., Vaughan, B. A., Michelson, P. F., Mitsuda, K., Lewin, W. 
H. G., van Paradijs, J., Penninx, W., & van der Klis, M. 1991, ApJ, 379, 295 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



- 16 - 



-400 -200 



200 400 



HJ_ 4 (*„ rb )! 2 



* rh =413.44 radians 



500 

Frequency (Hz) 



~K (*orb)! 2 




< r - r spin) / r orb 

-10 1 2 3 



- l J -3(*orb)! 2 


IJ_,( 


*„rb)! 2 


1-1,(4 


orb)!' 


U 3 (« 


orb)!' ' 










pin 














Jorb 

bins 


"orb 
bins 












* )l 2 


> 


* 

| 


|J 2 (4 


orb)! 2 




„A 






1 

IWmA 










-1 


o- 3 


-5x10 





5x1 0-* 


10 3 



Frequency-f spin (Hz) 



Fig. 1. — These two panels show the Fourier response of the fundamental spin harmonic of a phase modulated 
2 ms binary pulsar. The properties of the response are typical for circular orbits with 3> or (, = 2nx or tf sp i n 3> 1. The 
response was calculated for an 8 hour observation of a NS-WD binary with a 50min orbital period, projected semi- 
major axis x or b = 0.132 lt-s, orbital inclination i — 60°, and a companion mass of ~ 0.2MQ. The panel on the right 
shows the central portion of the full response on the left. The periodic nature of the sidebands with peak-to-peak 
spacing (in Fourier bins) of r or t = T b s /Porb is obvious. 
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Fig. 2. — Results of the measurement of the maximum and average values of J s {&orb)\ where s are all the integers 
from — > \&orb] ■ The Bessel functions correspond to the magnitude of the s th sideband pair of the Fourier response 
of a phase modulated sinusoid as a fraction of the unmodulated or fully coherent amplitude (see §2.1). Approximate 
asymptotic relations for large values of <J> ori) are given above the measured values. 
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Fig. 3. — Orbital versus spin periods for 54 currently known binary pulsars with P OT b < 10 d. The symbols denote 
the type of system (circles for NS-WD systems, squares for NS-NS systems, and stars for the three recently discovered 
accreting X-ray MSPs) as well as its location in the galaxy or a globular cluster (open or filled respectively). Lines 
of constant phase modulation amplitude $ or (, = 2-KX or bf sp in are overlaid assuming orbital inclinations of 60° , pulsar 
masses of 1.4 Mq, and companion masses of 0.1, 1, and 10 M for dotted, dashed, and solid lines respectively. 
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Fig. 4. — Histograms of spin period, orbital period, and phase modulation amplitude $ ori) for the 54 binary pulsar 
systems described in Fig. 3. 
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Fig. 5. — Two examples of key parts of a phase modulation sideband search. The left column shows the analysis of 
an 8 h observation of the globular cluster 47 Tucanae taken with the Parkes radio telescope of the binary millisecond 
pulsar 47TucJ. The right column shows analysis of simulated data containing a weaker millisecond pulsar. Plots a and 
e show portions of the full length power spectra P(r) centered on the fundamental spin harmonics r spin — f sp inT bs 
of the pulsars. The grey regions are the sections of the power spectra P(r) that were Fourier analyzed in order 
to create the power spectra Pzfr) shown in plots b and /. Plots c and g show histograms of the orbital phase as 
measured using the raw Fourier amplitudes in the original FFTs as described in §4.2.2. Once the three Keplerian 
elements for circular orbits were determined, a Fourier domain orbital template was calculated and correlated with 
the raw Fourier amplitudes (plots d and h). This process recovered all of the power from the pulsar's fundamental 
spin harmonic that had initially (plots a and e) been spread over many Fourier bins. The dotted lines in each plot 
show the known values of the independent variable. A more detailed description of each plot is given in §4.2. 
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Fig. 6. — The results of simulations showing the sensitivity of both acceleration and sideband searches for the 
fundamental harmonic (i.e. a sinusoid) of a 2 ms binary pulsar with a 0.1 Mq companion and orbital periods of 0.5, 
1, and 2h as a fraction of the optimal (i.e. coherent) sensitivity, which is proportional to T~^J 2 . The detection 
threshold used in the simulations was 8 a. Sideband searches are significantly more sensitive to binary pulsars than 
optimal duration acceleration searches when T b s > 2P or b- 
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Fig. 7. — Results from simulations similar to those shown in Figure 6 but for a 20 ms binary pulsar with a 1.0 Mg 
companion. 



